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A novel method is presented for the beam shaping of far field intensity distributions of coherently combined 
fiber arrays. The fibers are arranged uniformly on the perimeter of a circle, and the linearly polarized beams of 
equal shape are superimposed such that the far field pattern represents an effective radially polarized vector 
beam, or discrete cylindrical vector (DCV) beam. The DCV beam is produced by three or more beams that 
each individually have a varying polarization vector. The beams are appropriately distributed in the near 
field such that the far field intensity distribution has a central null. This result is in contrast to the situa- 
tion of parallel linearly polarized beams, where the intensity peaks on axis. © 2009 Optical Society of America 
OCIS codes: 140.3298, 140.3300, 060.3510. 
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The propagation of electromagnetic fields from multiple 
fibers can result in complex far-field intensity profiles 
that depend crucially on the individual near field polar- 
ization and phase. Coherently phased arrays have been 
used in defense and communications systems for many 
years, where the antenna ensemble forms a diffraction 
pattern that can be altered by changing the antenna 
spacing, amplitude, and relative phase relationships. In 
the optical regime, similar systems have recently been 
formed by actively or passively locking the phases of two 
or more identical optical beams [1-4]. While coherently 
combined fiber lasers are increasingly gaining acceptance 
as sources for high power applications, most approaches 
involve phasing a rectangular or hexagonal grid of fibers 
with collimatcd beams that arc linearly polarized along 
the same axis [5]. The resultant far field diffraction pat- 
tern therefore is peaked in the center. If the polarization 
is allowed to vary, as it does in radial vector beams, the 
diffraction pattern vanishes in the center [6], which can 
result in a myriad of device applications, discussed be- 
low. 

Several types of polarization states have been inves- 
tigated, including radial and azimuthal [6-18]. These 
beams are used in mitigating thermal effects in high 
power lasers [19-21], laser machining [22,23], and par- 
ticle acceleration interactions [24,25]. They can even be 
used to generate longitudinal electric fields when tightly 
focused [26,27], and although they are typically formed 
in free space laser cavities using a conical lens, or axicon, 
they can also be formed and guided in fibers [28-30]. 
Typical methods for creating radially polarized beams 
fall into two categories: the first involves an intracavity 
axicon in a laser resonator to generate the laser mode, 
the second begins with a single beam and rotates the po- 
larization of portions of the beam to create an inhomo- 
geneously (typically radially or azimuthally) polarized 
beam. In this Letter we demonstrate a novel approach 



that takes an array of Gaussian beams, each with the ap- 
propriately oriented linear polarization, and then super- 
impose them in the far-field, effectively creating a DCV 
beam. 

Our starting point in determining the propagation of 
polarized electromagnetic beams is the vector Hclmholtz 
equation, V x V x E(p, z) — fc 2 E(p, z) = 0. We use cylin- 
drical coordinates, so that p is the usual transverse coor- 
dinate, fco = io/c, and the usual sinusoidal time depen- 
dence has been factored out. Applying Lagrange's for- 
mula permits the wave equation to be written in terms 
of the vector Laplacian: V 2 E(p, z) + fc 2 E(p, z) = 0, 
which can then be reduced to the paraxial wave equation. 
The paraxial solutions arc then inserted into the familiar 
Fraunhofcr diffraction integral, which for propagation at 
sufficiently large z, yields the electric field, 

ikop 2 \ exp(ikoz) 
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where E° (r + Sj , 0) is the incident electric field in the 
plane of the fiber array. It is clear that longitudinal 
polarization is not considered here, consistent with the 
paraxial approximation. The far field intensity will be 
a complex diffraction pattern that depends on the in- 
dividual beam intensity profile in the near field as well 
its polarization state and optical phase. A diagram il- 
lustrating the geometry used is shown in Fig. 1, where 
an array of three circular holes of radius a are equally 
distributed on the circumference of a circle of radius R. 
The center of each hole is located at Sj = Rij, where 
Yj = (SccosOj + ysinOj) is the unit radial vector, and r 
is the coordinate relative to each center. The center of 
each fiber is separated by an angle 9j = 2ir(j — 1)/N. 
One can scale to any number of beams, constrained 
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We employ two different methods to create the DCV 
beams previously described. In each method, the experi- 
mental configuration utilizes collimated Gaussian beams 
that are in phase. The arrays are cylindrically symmet- 
ric both with respect to intensity and polarization. We 
focus the beams with a transform lens in order to sim- 
ulate the far field intensity at the lens focus. The focal 
spot is then imaged onto a camera by a microscope ob- 
jective in order to fill the CCD array. The laser source is 
a 100 mW continuous wave (CW) Nd:YAG operating at 
A = 1.064/im. 



Fig. 1. Example of the multiple fiber setup. A specific 
case of N = 3 is shown, where each beam is represented 
by a circle of radius a and arranged uniformly on a circle 
of radius R. Points within the zth hole are located by 
d = Si + r, with the vector r originating at the center of 
each hole. 
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only by the radius R. For identical Gaussian beams, 
the incident field is expressed as, E° = E°(r)£j, where 
J5°(r) = E e- r ^ w o . This permits Eq. (1) to be sepa- 
rated, 
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E(k X ,ky)=F( P , Z)Y,< 
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where F(p,z) = E ko/z rdr J (k rp/ z)e~ r l w °, and 
the prefactors that do not contribute to the inten- 
sity have been suppressed. It is clear from Eq. (2) 
that the far field transform preserves the radially 
symmetric polarization state of the system, as ex- 
pected. Note that in the limit wo/a 3> 1, and using 
the relation xJq(x) = [xJi(x)]' , we have !F(p,z) ~ 
Eo(a/ p)J±(koap/ z), which gives the expected Fraunhofcr 
diffraction pattern for a circular aperture of radius a. In 
the opposite limit, where the Gaussian profile is narrow 
enough so that the aperture geometry has little effect, 
we have !F{p,z) « E k WQ/ (2z) exp[—(k pwo/(2z)) 2 ]. 
For the multiple beam arrangements investigated, cor- 
responding to N = 3,4,6, we then can calculate the 
intensity, Tjv, explicitly: 
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where In = T (p, ^)/Cjv(</>; p, z), k x = (ko / z) p cos <j> and 
k y = (ko/z)psiii(f>. To rotate the array, one can perform 
a standard rotation r' = lZ(4>')r, so e.g., a ir/4 rotation 
would give, /C 4 (>') = 4(1 - cos(V2k x R) cos(y/2k y R)). 



Fig. 2. Diagram of the experimental setup for DCV beam 
generation. The oscillator output is split and phase mod- 
ulated (S/PM). The output of each fiber is collimated by 
means of a lens, LI, and then the polarization is rotated 
by a half- wave plate (HWP). The beam ensemble is then 
focused to a point which is imaged by a microscope ob- 
jective (MO) onto a charge-coupled-dcvicc (CCD) array. 

The first method (for the case of N = 4) involves 
expanding a beam from the Nd:YAG by means of a tele- 
scope beam expander to a diameter of about 1cm. The 
diffractive optical element then creates an 8 x 8 array of 
beams. Four of the beams are reflected and made nearly 
parallel by a segmented mirror, and then linearly polar- 
ized upon passing through a set of four half-wave plates. 
The phases of each beam are made identical by having 
them pass through articulated glass slides. 

The second method (for the case of N = 3 and 6), 
which is simpler for larger arrays, is shown in Fig. 2. In 
this method the Nd:YAG is coupled into a polarization 
maintaining fiber. This signal is then split into multi- 
ple copies by means of a lithium niobatc waveguide 8- 
way splitter and 8-channcl electro-optic modulator. Each 
path can then have a separate phase modulation to en- 
sure proper phasing of the beams. These signals are then 
propagated through fiber and coupled out and collimated 
by a 1 in. lens. The lens is slightly overfilled so the Gaus- 
sian outputs of the fibers are truncated at 1.1 x (the e _1 
radius). 

Experimental results for the radial vector beam gen- 
eration found good agreement with theoretical calcula- 
tions, as illustrated in Fig. 3. Clearly, the interference 
patterns and symmetry correlate well with the corre- 
sponding calculations in the three beam case. In particu- 
lar, for N = 4 (middle row), the central null surrounded 
by a checkerboard peak structure. This pattern differs 
from the linear polarization by a rotation of 7r/4, based 
on a simple phase argument. The N = 6 case is shown in 
the bottom row, where clearly the intensity vanishes at 
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Fig. 3. Far field intensity profiles of effectively radially 
polarized beams. The left panels correspond to the ex- 
perimental images and the right panels represent the 
theoretical results. We consider a 3 (top row) , 4 (middle 
row), and 6 (bottom row) beam arrangement. For N = 
we take wq = 32 pm and R = 5wq. For N — 4. wo = 42 
/im, and R = 2.7wo. For N = 6, we have wq = 43 pm 
and R = 3.8wq. Note that in general, a radially polarized 
beam undergoes a discontinuity at the origin, and thus 
the intensity must vanish there. 



the origin, followed by the formation of bright hexagonal 
rings. The three beam configuration is shown in the top 
row, where again there is satisfactory agreement between 
the calculated and measured results. Any observed dis- 
crepancies may be reduced with an appropriately incor- 
porated feedback system. A substantial fraction of the 
energy in the plots of Fig. 3 is contained within the first 
regions of high intensity peaks. We illustrate this for the 
N = 4 case where the net intensity within a circular re- 
gion of radius p, U(p), is calculated to be approximately 
U(p) w k pz/(2R)[k pR/z-Ji(2k pR/z)}. Inserting the 
appropriate parameters and normalizing to the intensity 
integrated over the entire image plane, yields 50% of the 
distribution is contained in the neighborhood of the first 
main peaks. 

In conclusion, it has been shown by employing two 
different configurations, and by carefully controlling the 
polarization, a central null can be created in the far field. 
In previous work, the central portion of a phased ar- 
ray of lasers typically has been a peaked function. Now 
for the first time to our knowledge, a central null has 
been formed. Similar results have been shown from sin- 
gle aperture lasers using an axicon to generate an annular 
mode inside a laser cavity. In the present case however, 



there is direct control over which mode will be generated 
from the same system. As for scaling, there is no funda- 
mental limitation on the number of beams used. Several 
methods are available to incorporate a greater number of 
beams, including concentric beam placement (following 
the same beam-combining prescription described in this 
Letter) that could scale in roughly the same pattern as a 
Bessel beam. It is also possible to use fiber lasers with ac- 
tive phasing, thus potentially scaling to many hundreds 
of beams [5]. Our system also has a practical advan- 
tage over typical high power Gaussian beam applications 
which have the drawback of having the beam concen- 
trated near the center, where the reflecting beam direc- 
tor (such as a telescope obscuration) resides. Our device 
on the other hand, has no power propagating along the 
center axis of the beam, and thus for high power ap- 
plications requiring a center obscuration telescope beam 
director, our proposed system offers new advances. 
We thank S. Feng for many useful discussions. 
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